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ABSTRACT 


The  introduction  of  complex  Navier -Stokes  equations  shows  that 
steady  axisymmetric  motions  of  viscous  incompressible  fluids  around 
conical  surfaces  can  be  expressed  in  terms  of  the  corresponding 
general  solution  of  the  Stokes  equations  of  slow  motions.  The  latter 
integration  is  accomplished  with  the  aid  of  slow-motion  eigenfunctions 
with  integral  eigenvalues  for  infinite  plates  and  semi-infinite  needles 
and  with  generally  complex  eigenvalues  for  cones  and  conical  corners. 
The  eigenvalues  and  eigenmotions  obtained  resemble  the  corresponding 
eigenvalues  and  eigenmotions  of  the  analogous  flows  past  dihedral 
angles.  In  particular,  the  existence  of  critical  and  branching  eigen¬ 
values  reveals  that  laminar  flows  past  conical  surfaces  depend  on  the 
cone  angle  in  a  nonanalytic  manner.  The  investigations  include  a  note 
on  diffusor  and  jet  flows. 
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1 .  Introduction 

Slow  motions  around  the  edge  of  a  dihedral  angle,  which  are  governed 
by  the  linear  Stokes  equations,  have  been  extensively  studied  in  [4,  6, 
and  10],  By  means  of  "complex  Navier-Stokes"  equations  introduced  in  [10] 
it  became  possible  to  reveal  four  fundamental  properties  of  viscous 
fluid  motions  which  are  important  for  the  general  theory  of  the  Navier- 
Stokes  equations: 

(I)  The  solutions  of  the  Navier-Stokes  equations  are  essentially 
representable  by  "slow-motion  eigenfunctions." 

(II)  Flow  patterns  of  laminar  motions  which  are  "attached"  or  "separated" 
at  a  sharp  or  blunt  edge  of  an  obstacle  can  be  studied  in  an  exact 
manner . 

(III)  In  contrast  to  inviscid  fluid  motions  viscous  fluid  flows  around 
obstacles  exhibit  "critical  body  parameters"  which  prohibit  analytic 
parameter  expansions  without  rigorous  justification. 

(IV)  The  Navier-Stokes  equations  yield  solutions  which  represent  "regular 
motions"  with  bounded  velocities  and  pressures, and  other  solutions 
which  describe  "singular  motions"  with  infinite  pressures  and  finite 
or  infinite  velocities. 

The  present  paper  confirms  the  properties  (I)  through  (IV)  by  investi¬ 
gating  the  analogous  problem  of  viscous  fluid  motions  around  a  conical 
surface.  For  the  sake  of  simplicity,  the  derivations  will  be  confined  to 
axisymmetric  flows  which  depend  on  two  variables.  The  general  problem 
considered  includes  as  a  special  case  the  symmetric  flow  past  a  semi¬ 
infinite  needle  which  has  been. briefly  studied  in  [9]. 
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2 .  Motions  Around  Conical  Surfaces 

Consider  the  steady  axisymmetric  motion  of  a  viscous  incompressible 
fluid  around  a  conical  surface  (Figs.  1  and  2),  the  solid  angle  of  which 
may  be  called  "cone  angle"  2a  (0  ^  2a  <  2tt).  The  complementary  angle 
2p (=  2tt  -  2a)  is  called  the  "corner  angle"  of  the  conical  surface.  For 
2a  =  0  the  conical  surface  degenerates  to  a  semi-infinite  thin  needle 
and  for  2a  =  rr  to  an  infinite  flat  plate.  For  0  <  2a  <  tt  the  conical 
surface  represents  a  cone  and  for  0  <  23  <  tt  a  conical  corner. 

Let  (r,9)  denote  a  polar  system  in  the  meridian  plane  of  the  axi¬ 
symmetric  motion  such  that  9  =  ±  tt  designates  the  solid  axis  of  the 
conical  surface  (Fig.  1).  Furthermore,  let  (u,v)  represent  the  corre¬ 
sponding  radial  and  tangential  velocity  components  of  the  motion  with 
the  variable  pressure  p.  If  r,u,v,  and  p  are  reduced  to  dimensionless 
variables  by  appropriate  flow  parameters  which  are  characteristic  for 
the  motions  considered  (see,  for  instance , [11  ]) ,  then  the  "real  Navier- 
Stokes  equations"  can  be  written  in  the  form: 


-  (ru)  +  -4-  u  +  (ufi  -  2v)  - 

r  rr  r3  09  r2  9 


—  (vQ  +  u)  =  f  p„  +  uu 
r 


.2  '*9  '  5  1  r  r 


+  J  (u0  -  v)  , 


r  (rv)rr  +  ^  v99  +  F  (v  cot0)0  +  ^  u9  =  £  p0  +  uvr 


+  7  (v0  +  u)  , 


(urs  sin9)  +  (vr  sin0)Q  =  0. 

V  u 


(1) 

(2) 

(3) 
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The  latter,  continuity  equation,  leads  to  the  introduction  of  the 
stream  function  \|r  which  is  defined  by 


u  = 


r2sin0 


v  = 


r  sin9  Tr 


(4) 


As  was  shown  in  [10]  it  is  advantageous  to  introduce  the  concept  of 
"complex  laminar  flows"  which  are  described  by  the  "complex  velocity" 
(U,V),  the  "complex  pressure"  P,  and  the  "complex  stream  function" 

¥ ,  the  real  parts  of  which  represent  "real  laminar  motions": 


u  -3L 


u 


v  =  2- 


P  = 


P  +  P 


Y  +  Y 


(5) 


It  is  easily  verified  that  complex  laminar  motions  are  governed  by  the 
following  "complex  Navier-Stokes  equations" 


(rU)  +  \  U  +  (U  -2V)  ■  4  (V  +D)  *  i  P 

rr  r2  06  r“  0  r3  q  ~  x 


‘  4 


(U+U)Ur  +  U(Ur+Ur)  +  (Uq-V)  +  2  (Uq+Uq-V-V) 


<rV>  +  -V  Vno  +  h  (y  cote),  +  2_  u  p 


rr  r2  99  r2 


e  t.2  e  2r  e 


+  z|^(U+U)Vr  +  U(Vr-fVr)  +  — ±-V  (Ve+U)  +  ^  (Vq+Vq+U+U)  j  , 


(Ur3sin0)  +  (Vr  sin6)  =  0, 
r  0 


(6) 

(7) 

(8) 


which  differ  from  the  real  Navier-Stokes  equations  only  in  their  non¬ 
linear  inertial  terms. 
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Axisymmetric  complex  motions  around  conical  surfaces  are  defined 
as  those  solutions  of  the  complex  Navier-Stokes  equations  which  satisfy 
the  symmetry  condition 


U(r,-0)  =  +  U (r  j-l-0)  ,  V (r , -6)  =  -  V(r,+0)  (9) 

and  the  nonslip  condition 

r^O,  6  =  ±  3  :  U  =  0,  V  =  0  .  (10) 

When  a  complex  solution  is  found,  the  corresponding  real  motion  is 
obtained  by  deleting  the  imaginary  parts  which  have  no  physical  meaning. 
With  the  complex  stream  function  ¥  defined  by 


U  - - v 


V 


¥ 


r2sin0  6  *  r  sinG  r 

the  complex  Navier-Stokes  equations  (6)  and  (7)  can  be  combined  to  the 
single  equation 


(11) 


E4Y  =  MOO 


(12) 


for  the  stream  function  f.  The  symbol  E3  denotes  the  linear  differential 
operator 


E3Y  =  Y  +  —  Y 


rr 


89 


cot_9 


(13) 


and  M'  the  nonlinear  differential  operator 


M(f)  = 


siM‘U+7)  f-Jfl  )  -  cm)  /,  \ 

L  r\r2sin20/G  e\r2sin30/r 


(14) 
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For  real  stream  functions  i|r  the  latter  nonlinear  operator  M  reduces  to 
the  well-known  form 


M(i)r)  =  sin8 


E3t  \  _  t  /  E2fr 

2 sin3 8  /  8  6  \  r2sin2  8 


),] 


Since  the  Navier-Stokes  equations  are  linear  in  their  derivatives 
of  highest  order,  they  may  be  considered  as  "almost  linear"  partial 
differential  equations.  Hence,  as  was  pointed  out  in  [10]  exact 
solutions  can  be  constructed  in  an  absolutely  Linear  manner  by  solving 
at  first  the  linear  homogeneous  Stokes  equations  of  slow  motions  (16), 
(17),  and  (18)  which  neglect  all  inertial  terms  of  the  Navier-Stokes 
equations.  Subsequently,  a  solution  of  the  complete  equations  can  be 
found  by  integrating  successively  inhomogeneous  linear  equations  which 
result  from  the  basic  slow-motion  solution.  A  detailed  description 
of  the  latter  integration  step  is  presented  in  [10],  In  order  to  avoid 
a  formal  repetition,  only  the  general  slow-motion  integral  may  be  con¬ 
sidered  further  because  it  retains  the  exact  properties  of  the  flow 
around  the  vertex  of  the  conical  surface  (see  [10]). 

3.  Separable  Slow-motion  Integrals 

According  to  the  previous  section  the  (complex)  Stokes  equations 
of  slow  motions 


r  +  js  uee  +  ^  <V2V)  ’  7  <V“>  '  *5r  " 


-  (rV) 
r  v  'rr 


+  i-  V 


+  ~  (V  cote)n  +  2-  U 


-L-  P 


28  2r  0  ’ 


5 


NWL  REPORT  NO.  1922 


can  be  combined  to  the  differential  equation 


E  Y  -  Yrrrr  +  ^2  ^rr00 


-  cot0Yrr0)  -  jz  (Yr0e  -  cot0Yr0) 


+  —  v 
r4  ”0000 


-  l£?fs  +  ir  (8  +  3cotse» 


00 


-  (3  +  cot2  0)Y0  =  0 


(18) 


for  the  stream  function  Y.  Although  the  differential  equation  (18)  is 
not  separable  in  the  usual  sense,  "separable  solutions" 


Y  *  G (r)H(0)  (19) 
can  still  be  found  in  the  same  generalized  manner  shown  in  [10],  A  sub¬ 
stitution  of  Eq.  (19)  into  Eq.  (18)  leads  to  the  equation 


H/y//  -  2cot0Hw  +  [8  +  3cot2 0  +  |  (r2G"  -  2rG/)]H// 

0  n  mr 

-  [3cot0(3  +  cot2 0)  +  cote  £  (r2G//  -  2rG')]H'  +  r4  H  =  0,  (20) 

which  splits  into  the  three  ordinary  differential  equations 
r2G"  -  2rG'  -  u^G  =  0  ,  (21) 

r4G////  -  0UgG  ®  0  ,  (22) 

U""  -  2cot0H///  +  [8  +  3cot2 0  +  2^  ]H" 

-  [3(3  +  cot2 0)  +  2ui1]cot0H/  +  (ju2H  =  0  ,  (23) 

where  a)!  and  0)3  are  arbitrary  complex  constants. 


The  first  two  equations  are  of  Eulerian  type  and  can  be  integrated 
by  functions  of  the  form 

A+2 


G  =  r 


(24) 
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The  corresponding  characteristic  equations 

=  (X  -  1)(X  +  2),  uig  =  (X  -  l)X(X  +  1)(X  +  2)  (25) 

yield  the  relationship 

aus  =  X(X  +  Da^  (26) 

between  and  xs ,  so  that  they  both  may  be  replaced  by  the  exponent  X. 
Since  the  velocity  must  be  zero  at  the  vertex  of  the  conical  surface, 
the  real  part  of  X  must  fulfill  the  "weak  regularity  condition"  (see 
Eq.  (11)) 

X  +  X  >  0  .  (27) 

A  stronger  regularity  condition  will  be  required  later  on  the  basis  of 
a  bounded  pressure  (see  [10]). 

The  remaining  fourth  order  differential  equation  (see  [1]) 

H""  -  2cot0H"/  +  (3cot2 0  +  2X2  +  2X  +  4)H"  -  (3cot20  +  2X2  +  ZX  +  5)cot0H' 
+  (X  -  1)X(X  +  1)(X  +  2)H  =  0  (28) 

can  be  factorized  in  a  commutative  way: 

[3$  '  cot0  de+  x(x  '  l)][de?‘  cote  fi+  a  +  l)a  +  2,]H  ‘  0>  (29) 

Hence,  one  has  the  two  second  order  differential  equations  (X  ^  -  -g-) 


H"  -  cotOH"  +  (X  +  1 ) (X  +  2 )H  =  0,  (30) 

U"  -  cot0H/  +  (X  -  1)XH  =  0,  (31) 

which  assume  the  simple  forms 

(1  -  t2 )H  +  (X  +  1) (X  +  2 )H  =  0,  (32) 

(1  -  1 2)ii  +  (X  -  1)XH  =  0,  (33) 


where  the  dots  denote  differentiation  with  respect  to  t  =  cos0.  The 
product  form  (29)  of  the  differential  equation  (28)  is  not  surprising 
because  the  differential  equation  (30)  yields  the  well-known  potential 
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solutions  of  the  same  problem  which  represent  simultaneously  slow 
motions . 

It  is  easily  verified  that  the  differential  equations  (32)  and  (33) 
are  integrable  with  the  aid  of  the  first  order  associated  Legendre 
functions  of  first  and  second  kind 


P*(t)  =  / 1  -  t2  Pv(t),  Q^(t)  =  / 1  -  t2  Qv(t)  ,  (34) 

where  the  corresponding  degrees  v  =  X  +  1  and  v=X-1^0of  the 
Legendre  functions  Pv(t)  and  Qv(t)  may  be  real  or  complex.  Thus,  one 
arrives  at  the  following  set  of  separable  slow-motion  integrals  of  the 
linear  differential  equation  (18): 

Y  =  rX+2sin0[A  .  P*  (cos0)  +  B  P*  (cos0)  +  C  (cos0)  +  D  (cos0)] 
A  A  A-rl  A  A  "I  A  A"*"!  A  A—l 

(35) 


for  X  +  X  >  0  and  X  ^  1,  and 


=  r3sin9[A1Pg  (cos0)  +  Bx  tan  |  +  CxQ^(cos9)  +  Da  ] 


(36) 


for  X  ■  1.  The  coefficients  A  ,  B  ,  C.  ,  and  D  are  arbitrary  complex 

AAA  A 

numbers.  For  vanishing  coefficients  B^  and  the  slow  motion  solutions 
(35)  and  (36)  represent  potential  flows. 

4.  Symmetric  Slow-motion  Eigenfunctions 

The  separable  slow-motion  integrals  (35)  and  (36)  represent  axi- 
symmetric  slow  motions  around  a  conical  surface  if  and  only  if  the  sym¬ 
metry  condition  (9)  and  the  nonslip  condition  (10)  are  satisfied.  Un¬ 
fortunately,  the  Legendre  functions  of  first  and  second  kind  Pv(t)  and 
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Q^(t)  are  generally  singular  at  t  =  ±  1  (9  =  0  and  9  =  tt)  .  Consequently, 
it  is  more  advantageous  to  represent  the  slow-motion  integrals  (35)  by 
the  associated  Legendre  functions 

^(t>  =  t2iv(t>,  ai(t)  =  y  i  -  t2uv(t>  07) 

in  which  the  Legendre  function  of  first  kind  $v(t)  is  regular  at 
t  =  +  1  (9=0)  for  all  degrees  v. 

Since  the  Legendre  function  of  second  kind  fc^(t)  is  singular  at 
t  =  +  1  (9=0)  for  all  degrees  v,  the  " ax i symmetric  slow-motion  eigen¬ 
functions"  of  conical  surfaces  are  of  the  form 

=  rX+2[AxHx+1(cos9)  +  B^H^  (cos9)  ]  (38) 

for  \  +  \  >  0  and  \  1 ,  and 


\  =  r^Aj.Hg  (cos9)  +  Bx  (1  -  cos9)]  (39) 

for  \  -  1.  While  the  coefficients  A.,  and  B,  are  complex  numbers  for 

A.  A. 

complex  exponents  they  can  be  confined  to  real  values  for  real  ex¬ 
ponents.  The  term  Hv(t)  is  used  as  an  abbreviation  for  the  function 

Hv(t)  =  /rr?*;(t)  =  a  -  ts)  *Bv(t)  (v  ^  o).  (40) 

One  derives  from  Eqs.  (32)  or  (33)  the  following  recurrence  formula 


h  _  n(n  -  1)  -  v (v  +  1) 

v’n+1  "  2n(n  +  1) 

for  the  coefficients  of  the  expansion 


hv>n  (n  1>2>3’  •••) 


(41) 


H  (t)  =  2  h  (1  -  t)n. 

v  n  =  1  v»n 


(42) 
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After  combining  Eqs.  (11),  (16),  (17),  (38)  and  (39)  one  arrives 
at  the  velocity  (U^,  V^)  and  the  pressure  P^  (neglecting  an  arbitrary 


pressure  constant)  of  the  slow  motions  (38)  and  (39): 

=  r^[A^H^+1(cos9)  +  B^_1  (cos0)  ]  ,  (43) 

v  =  A-±_2  rX[A  H  ..  (cose)  +  B  H  (cos 6)]  ,  (44) 

A  sin0  A.  A  •  1  A  A“1 

P  =  4  rX_1H  (cose)  ,  (45) 

A  A  -  1  A~1 

for  A  +  A  >  0  and  A  /  1,  and 

Ux  =  r[ALE3(cos0)  -  Bx  ]  ,  (46) 

Vi  =  iife  ^AM0080)  +  V1  -  cos0)]  ,  (47) 

Px  =  -  12B1  log (r  cos2  |  )  (48) 

for  A  =  1. 


The  behavior  of  the  pressure  functions  (45)  and  (48)  shows  that  the 
concept  of  "regular  flows"  and  "singular  flows",  which  was  introduced  in 
[10]  according  as  the  pressure  of  the  motion  is  bounded  or  unbounded, 
applies  also  to  motions  around  conical  surfaces.  In  fact,  one  has  the 
following  theorem  for  conical  surfaces  which  is  equivalent  to  the 
corresponding  theorem  for  dihedral  angles. 

THEOREM  1:  Eigenfunctions  of  regular  axisymmetric  motions  around 

conical  surfaces  contain  only  "regular  eigenvalues"  A  which  fulfill 
the  "strong  regularity  condition" 

X  >  1  .  (49) 

2 
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The  term  "strong"  regularity  condition  is  justified  in  contrast  to 
the  notion  "weak"  regularity  condition  which  was  used  for  Eq.  (27). 

In  opposition  to  other  opinions  (see,  for  instance,  [5])  it  was 
pointed  out  in  [7,  9, and  10]  that  only  regular  motions  appear  to  be 
physically  realistic  as  they  exhibit  the  features  of  actually  observed 
flows.  Furthermore,  in  a  paper  in  preparation  it  will  be  demonstrated 
by  means  of  compressible  fluids  that  the  acceptance  of  singular  motions 
leads  to  such  objectionable  consequences  as,  for  instance,  positive  or 
negative  infinite  densities. 

The  nonslip  condition  (10)  and  Eqs .  (43),  (44),  (46),  and  (47) 
lead  to  the  homogeneous  algebraic  equations  which  determine  the  coef¬ 
ficients  A  and  B  of  the  slow-motion  eigenfunctions  (38)  and  (39): 

X  X 


Wi(cose)  +  BxVi<cose) ' 0 

A  H  (cos0)  +  B  H  (cosS)  =  0 
X  X+l  X  X-l 

for  X  +  X  >  0  and  X  ^  1 ,  and 

A1H2(cos0)  +  Bj  (1  -  cos3)  =  0 
A1H2(cos0)  -  Bx  =0 

for  X  =  1 . 


(50) 


(51) 


With  Hg (t)  =  (1  -  t2)t  the  algebraic  system  (51)  yields  nontrivial 
solutions  and  Bx  if  and  only  if  the  corner  angle  20  =  20*  satisfies 
the  "cosine  condition" 

* 

cos30*  =  cos2  j  ,  (52) 

which  leads  to  0*  =  120°. 
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THEOREM  2:  The  "singular  eigenvalue"  X  =  1  of  ax i symmetric  motions 
around  conical  surfaces  exists  only  for  the  "limiting  cone"  with 
the  cone  angle  la*  -  120°. 

The  critical  angle  of  the  limiting  cone  may  be  compared  with  the 
critical  angle  2of*  «  102°  of  the  limiting  wedge  which  was  found  in 
[10].  Because  the  limiting  wedge  angle  is  smaller  than  the  limiting 
cone  angle,  a  cone  appears  to  be  more  critical  than  a  wedge.  This 
is  physically  plausible  since  the  smaller  "volume"  of  a  cone  causes 
less  displacement  effects  in  a  motion  than  the  larger  "volume"  of  a 
wedge.  As  was  pointed  out  in  [9 and  10],  in  a  fluid  motion  the  dis¬ 
placement  effects  caused  by  the  "volumes"  of  a  blunt  wedge  or  cone 
(2q?“  <  2q?)  dominate  the  friction  effects  which  are  caused  by  the 
"surfaces"  of  the  wedge  or  cone.  At  a  sharp  wedge  or  cone  (2q-  <  la*) 
one  observes  the  opposite:  the  "surface"  friction  effects  exceed 
the  "volume"  displacement  effects. 

The  real  slow-motion  eigenfunction  for  the  limiting  cone 
(2q-  =  2<y*  =  120°)  with  the  singular  eigenvalue  1  =  1  is  (a,  =  real 
parameter) 

=  a1r3[sin9  sin28  +  sin3  ^  ]  .  (53) 

The  algebraic  system  (50)  yields  nontrivial  solutions  (A^,  B^)  if 
and  only  if  the  Wronskian  determinant 

W(t)  =  H^OOH^Ct)  -  H^+1  (t)H^_x  (t)  (54) 
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of  the  two  functions  H  (t)  and  H  (t)  vanishes  at  t  =  cosp.  Thus, 

A+l  A”1 

the  condition 

W(cosP)  =  (cos3)H^_j  (cosP)  -  (cosP)H^_-^  (cosP)  =  0  (55) 

determines  the  "eigenvalues"  X  ^  1  of  axisymmetric  slow  motions  around 
a  conical  surface  with  the  corner  angle  2p. 

After  differentiating  Eq.  (54)  and  applying  Eqs.  (32)  and  (33)  one 

has 

"(t)  -  Vh  -  ViVi  * 2  fAr  Vx(t)Vi(t)-  <56> 


With  W(l)  =  0  the  eigenvalue  condition  (55)  becomes  equivalent  to  the 
"orthogonality  relation" 

.  l 


W(cosp) 


--2(2X  +  1)J  Hx+1(t)Hx.l(t)  -0 


(57) 


cosp 


Thus,  the  relations  (40)  lead  to  the  transformed  orthogonality  condition 


W(cosp) 


-2<2X  +  I)f  (t)  (t)dt  =  0 


(58) 


COSP 


for  the  associated  Legendre  functions  !B1  (t)  and  (t) . 

X+l  X-l 

THEOREM  3:  The  eigenvalues  X  ^  1  of  axisymmetric  slow  motions  around  a 

conical  surface  with  the  corner  angle  P  are  the  arithmetic  means  of 

the  degrees  of  those  pairs  of  first  order  associated  Legendre 

functions  of  first  kind  3J1  (t)  and  3J1  ,  (t)  which  are  orthogonal 

X+l  X~1 

over  the  range  cosp  <  t  <  1 . 
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It  is  well  known  (see  Eq.  (41))  that  the  Legendre  functions  of  first 
kind^v(t)  (which  are  regular  at  t  =  +  1  (6  =  0)  by  definition)  assume 
regular  properties  at  t  =  -  1  (0  =  rr)  and  degenerate  to  the  Legendre 
polynomials  Pv(t),  if  and  only  if  the  degree  v  is  an  integer.  Further¬ 
more,  the  first  order  associated  Legendre  functions  of  integral  degrees 
(X  =  n)$^+1(t)  =  In+l^)  and  ^n-l^  =  Pn-1  ^  are  ort^°§onal  over  the 
range  (-  1  ^  t  <:  +  1).  Since  their  product  is  an  even  function  of  t 
for  all  n,  they  are  also  orthogonal  over  the  range  (0  t  <  +  1). 

Hence,  the  integers  X  =  n  =  2,3,4,  ....  are  eigenvalues  of  axisymmetric 
slow  motions  normal  to  an  infinite  plate  (2a  =  20  =  tt)  and  parallel  to 
a  semi-infinite  thin  needle  (0  =  tt)  .  Vice  versa,  it  is  not  difficult 
to  prove  that  the  integers  X  =  n  =  2,3,4,  ....  are  the  only  eigenvalues 
of  axisymmetric  slow  motions  past  an  infinite  plate  and  around  a  semi¬ 
infinite  needle.  Indeed,  for  the  infinite  plate  problem  the  analyticity 
of  any  slow-motion  eigenfunction  (38)  at  the  origin  r  =  0  follows  from 
the  analyticity  of  the  boundary  conditions  (see  [3  and  5]).  Hence,  the 

Legendre  functions  (t)  and  33  (t)  must  degenerate  to  the  Legendre 

a“Kl  a."1 

polynomials  pn+i (fc)  and  pn_1(t)  of  the  integral  degrees  X  =  n  =  2, 

3,4,  ...  .  The  same  follows  directly  for  the  semi-infinite  needle 

problem  as  the  eigenfunctions  (38)  must  all  vanish  at  the  needle 

0  =  0  =  TT. 

With  the  eigenvalues  of  axisymmetric  slow  motions  past  an  infinite 
plate  X  =  n  =  2,3,4,  ....  one  obtains  the  corresponding  real  eigen¬ 
functions  (ajj  =  real  parameters): 
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.  n+2 

* .  =  a  r  £ 


r 


(60) 


The  eigenvalues  of  axisymmetric  slow  motions  around  a  semi-infinite 
thin  needle  X  =  n  =  2,3,4,  ....  lead  to  the  corresponding  real  eigen¬ 
functions  (an  =  real  parameters): 


i|rn  =  anrn+^sin0 


pn+l<cos0)  " 


Pn+1<~  » 

pi-i<- 


Pn-1  (cose) 


(61) 


This  result  confirms  the  integral  which  was  found  in  [9]  for  the  same 
needle  problem  by  Taylor  series  expansions.  In  particular,  one  concludes 
from 

un  " 

(62) 

+  anrns in2 6  Pn+1(cose)  -  vg±1|_  ^  ^(cose) 

n-1 '  ' 


P 


2anrncos0  Pn+1(cos0)  -  - 


n+1 


(-  1) 


Vl<-  P) 


Pn-1 (cos0) 


'] 


that  all  radial  velocity  components  un  vanish  at  0=0.  Hence,  as  was 
pointed  out  in  [9]  any  viscous  fluid  passes  along  an  infinitely  thin  needle 
(which  has  neither  volume  causing  displacement  nor  surface  causing 
friction)  in  an  undisturbed  manner  without  producing  any  shear  stress 
at  the  needle. 

Since  the  first  order  associated  Legendre  functions  of  first  kind 
and  integral  degree  P^+^  (O  and  P^_^  (t)  are  only  orthogonal  over  the 
ranges  (0  ^  t  s  1)  and  (-  1  <  t  ^  +  1),  one  has  the  following  theorem. 
THEOREM  4:  Regular  integral  eigenvalues  X(^  1)  of  axisymmetric  slow 
motions  around  a  conical  surface  exist  if  and  only  if  the  conical 
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surface  degenerates  to  an  infinite  plate  (2q?  =  20  =  tt)  or  a  semi¬ 
infinite  needle  ( 2a  =  0)  in  which  cases  all  eigenvalues  coincide 

with  the  integers  X  =  n  =  2,3,4,  . 

For  axisymmetric  slow  motions  around  cones  (0  <  2a  <  tt)  and  into  conical 
corners  (0  <  20  <  tt)  the  leading  eigenvalues  X  f  1  have  been  computed  by 
solving  Eq.  (55).  This  has  been  accomplished  by  truncating  the  series 
expansion  of  W(cos0)  after  some  appropriate  term  of  N-th  order  in 
(1  -  cos0).  With  the  recurrence  formula  (41)  the  remaining  polynomial 
.  of  order  (2N  -  3)  in  X 

W(cos0)  «  -  S  1  (1  -  cos0)n(n  -  2j  +  l)h.  ,h  =  0  (63) 

n  —2  i=l  X+1,J  X-i,n-J+i 

retained  sufficient  accuracy  for  the  leading  eigenvalues  of  interest. 
Numerical  results  are  presented  in  Tables  1  and  2  and  in  Fig.  3. 

The  numerical  results  confirm  theorem  4  and  show  that  almost  all 
of  the  slow-motion  eigenfunctions  of  cones  (0  <  2oi  <  tt)  and  of  conical 
corners  (0  <  20  <  tt)  are  characterized  by  complex  eigenvalues.  A 
comparison  with  the  results  presented  in  [10]  reveals  that  the  eigen¬ 
values  of  symmetric  slow  motions  past  dihedral  angles  and  around  conical 
surfaces  display  essentially  the  same  involved  structure.  In  particular, 
the  existence  of  branch  points  shows  that  the  eigenvalues  of  slow  motions 
around  dihedral  angles  and  conical  surfaces  depend  on  the  corresponding 
wedge  or  cone  angles  in  a  nonanalytic  manner. 

It  is  easily  verified  that  the  real  eigenfunctions  of  slow  motions 
past  conical  surfaces,  which  are  characterized  by  complex  eigenvalues, 
describe  eddy-type  flows  of  the  same  sort  as  the  analogous  eigenfunctions 
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of  slow  motions  around  dihedral  angles  (see  [4,  6,  and  10]).  However, 
the  general  real  solution,which  is  the  real  part  of  the  linear  com¬ 
bination  (see  [2,  3,  5,  and  11]) 

Y  =  ZY  (64) 

A. 

where  the  sum  must  be  extended  over  all  eigenvalues  X  with  non¬ 
negative  imaginary  parts,  includes  also  the  eddy-free  motions  of 
practical  interest. 

5.  Note  on  Diffusor  and  Jet  Flows 

Since  the  product  (29)  is  commutative,  the  complex  exponents  X  are 
symmetric  to  the  point  X  =  -  1/2.  Hence  if  X  is  an  eigenvalue  of  a 
conical  surface,  which  satisfies  the  strong  regularity  condition 
X  +  X  >  2  for  the  vertex  r  =  0,  then  p,  =  -  (X  +  1)  is  an  eigenvalue 
of  the  same  surface  which  fulfills  the  strong  regularity  condition 
ji  +  |i  <  -  4  at  infinity  (r  =  °°) ,  and  vice  versa.  The  corresponding 

eigenfunction  Y  (r,0)  =  r3Y  (  —  ,0)  represents  a  "regular"  slow 

M-  X  r 

motion  which  has  neither  sources  nor  sinks  at  infinity. 

If  one  applies  the  same  "inversion  procedure"  (see  [10])  to  a  slow- 
motion  eigenfunction  of  a  conical  surface  which  is  weakly  singular  at 
the  vertex  r  =  0  (0  <  X  +  X  ^  2),  then  one  arrives  at  a  slow-motion 
eigenfunction  which  describes  a  source  or  sink  flow  at  infinity  (r  =  °°) . 
This  reveals  the  fact  that  weakly  singular  motions  around  sharp  cones 
with  angles  2a  ^  2a  =  120  (see  Theorem  2)  are  truly  "singular"  as 
they  exhibit  a  source  or  sink  of  infinitesimal  strength  at  the  vertex 
of  the  cone  (see  [9  and  10]). 
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The  reversible  inversion  property  -of  regular  motions  past  conical 

surfaces  is  only  unilateral  for  weakly  singular  motions.  Indeed,  while 
the  singular  eigenvalue  X  =  1  exists  only  for  the  limiting  cone  with  the 
angle  2<y*  ■  120°,  the  eigenvalue  (j,  =  -  (\  +  l)  =  -2  exists  for  all 
conical  surfaces  (0  £  2<y  <  2m)  .  The  corresponding  slow-motion  eigen¬ 
function  is  (a_2  a  real  parameter) 

\|f  ■  2a.  (cos0  -  l)[(cos0  +  l)cos0  -  (3co£»2g  -  1)],  (65) 

— 2  —2 

which  differs  from  Ackerberg's  solution  presented  in  [1]  by  an  additive 
constant  that  adjusts  the  stream  function  f_2(r,0)  to  the  conventional 
condition  t„2(r,0)  =  0.  It  may  be  noted  that  the  slow-motion  eigen¬ 
function  (65)  excludes  the  possibility  of  Ackerberg's  non-existence  paradox 
for  the  infinite  flat  plate  (20*  =  20  =  tt)  .  In  fact,  the  general  solution 
Y  of  the  complex  Navier-Stokes  equation  (12),  which  is  bounded  at  r  =  «=, 
is  of  the  form  (see  [10]) 

Y  =  Yg(r,8)  +  YM(r,0)  (66) 

where  Yg  is  the  general  solution  of  the  analogous  slow-motion  problem 
containing  the  r-independent  eigenfunction  (65)  and  all  eigenfunctions 
that  are  regular  at  r  =  ®.  Accordingly,  Ackerberg's  solution  (see  [1]) 
is  in  error  because  the  function  Y^(r,8),  which  adjusts  Yg  to  a  solution 
of  the  full  Navier-Stokes  equation  (12),  cannot  degenerate  to  a  finite 
sum,  unless  it  vanishes  identically  together  with  Yg  =  0  (see  [10]). 
Furthermore,  Ackerberg  based  his  solution  only  on  the  slow-motion 
solution  (65)  and  failed  to  observe  the  existence  of  an  infinite  number 
of  slow-motion  eigenfunctions  which  are  all  regular  at  r  =  ®.  As  was 
pointed  out  in  [10]  for  0  <  2a  <  tt  and  0  <  20  <  tt  the  function  Y^ 
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contains  only  terms  of  the  separable  form  (38)  because  the  eigenvalues 
(j,  =  -  (A  +  1)  constitute  an  arithmetic  progression  only  for  2a  =  0  and 
2cy  =  2(3  =  tt  (see  Theorem  4).  In  the  latter  case  the  fundamental  system 
of  separable  slow-motion  solutions  (38)  must  be  extended  by  the  so- 

called  "associated"  separable  slow-motion  solutions  of  order  m 
m 

¥a  =  sa5  ya  (m  =  1»  2’  3j  •••  )j  (67) 

which  were  introduced  in  [10]. 

For  example,  the  associated  separable  solution  of  first  order  is 
of  the  form  (see  Eq.  (38)): 

yJ  =  rX+2log  r[A^H  (cos 9)  +  B*H  (cos0)] 

A  A  A+l  A  A“I 


(68) 


+  r^CA1  H  (cosG)  +  B1  —  H  (cos  9)]. 

A  9A  A+l  X  SA  A-l 


In  general,  the  solutions  (67)  involve  the  functions  (log  r)m 

(m  =  0,  1,  2,  ...  )  and  the  "associated  functions" 
m 

Hv  "  Bv®  Hv(cos0)  (m  =  1,  2,  3,  ...  )  .  (69) 

The  latter  functions  satisfy  the  recurrence  relation  (m  =  1 ,  2,  3,  ...  ) 

_  ..m  ,  .  m  ,  .  m-1  .  .  m-2 

(l-ts)Hv  +  v(v+l)Hv  +  m(2v+l)Hv  +  m(m-l)Hv  =  0  .  (70) 

For  instance,  the  first  order  associated  function  H^(t)  is  a  particular 
solution  of  the  inhomogeneous  differential  equation 

(l-t2)ii^+  v(v+l)H*  -  -  (2v+l)Hv  ,  (71) 

which  is  regular  at  t  =  1.  Since  the  inhomogeneous  portion  of  Eq.  (71) 
is  a  homogeneous  solution  of  the  reduced  equation,  no  solution  H'*' 
degenerates  to  a  polynomial  even  if  the  exponent  v  is  an  integer. 
Nevertheless,  an  exact  solution  can  be  found  in  the  form 
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H1  =  S1(t)  + 

v  v 


V 

Ivl 


Hv(t)  log (1+t) 


(v  jt  0,  -  1), 


(72) 


in  which  reduces  to  a  polynomial  of  order  n  =(|v  +  1/2 |  -  1/2) if  the 
exponent  v  is  an  integer  v  =  l,±2,±3,  ...  .  After  substituting 
Eq.  (72)  into  Eq.  (71)  and  applying  Eq.  (40) one  arrives  at  the  dif¬ 
ferential  equation  (v  ^  0,  -  1) 


(i-t2)s'i  +  v(v+i)s?: 


=  2  (1-t) 

Ivl 


.(v+1^  -  (|vffe|+fc)t* 


gl  j 


(73) 


in  which  for  v  =  l,±2,±3,  ...  the  inhomogeneous  portion  reduces 
to  a  polynomial  of  order  n  =  (|v  +  1/2  |  -  1/2)  which  is  lower  than  the 
order  of  the  homogeneous  polynomial  solution  Hv(t).  Hence,  the  dif¬ 
ferential  equation  (73)  yields,  indeed,  a  polynomial  solution  of  the 
order  stated  above. 

According  to  these  derivations  the  general  solution  of  the  full 
Navier-Stokes  equations,  which  describes  laminar  motions  along  an 
infinite  flat  plate  (,2a  =  2|3  =  tt)  with  an  orifice  at  the  origin 
r  =  0  and  with  a  source  or  sink  of  finite  strength  at  infinity 
(r  =  “),  is  of  the  form  (see  [10]): 

t  "  s 0_  (r,e),  0  =  r2"n  n?  F.(e)(log  r)j  . 

n=z  11  11  j  =o  J 


(74) 


After  some  algebra  one  arrives  at  the  following  two  leading  terms  of  the 
asymptotic  expansion  (74)  (a_^  an<^  a-3  are  real  parameters): 

$-2  t.2  =  2a-2^t3  “  D  (t  =  cos0),  (75) 


0-3  =  V  (l-t)[(l+t)(a_3+9a32  log  ^  )  -  a22  (3+2t+2t2)]  . 


(76) 
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It  may  be  noted  that  the  analogous  needle  ( 2a  =  0)  and  limiting  cone 
(2a  ■  2a  ■  120°)  problems,  which  also  seem  to  have  caused  certain  dif¬ 
ficulties  in  Ackerberg's  investigations  [1],  are  solvable  without  the 
knowledge  of  the  associated  separable  slow-motion  solutions  (67). 

Since  the  slow-motion  solution  (65)  is  the  leading  term  in  any 
asymptotic  expansion  of  a  laminar  flow  along  a  conical  surface  with  an 
orifice  at  the  vertex  r  =  0  and  a  source  and  sink  distribution  of  finite 
strength  at  r  =  it  is  useful  to  investigate  the  corresponding  flow 
pattern.  As  the  velocity  components 

u_2  =  6a_2  pg-  (cos2 0  -  cos2P),  v_2  =  0  (77) 

show,  one  has  a  strictly  radial  dif fusor-type  flow  for  all  corner  angles 
0  <  2(3  <:  tt  (see  Fig.  4a).  However,  for  all  cone  angles  0  <  2a  <  n  the 
motion  assumes  the  character  of  a  jet-type  flow  as  shown  in  Fig.  4b.  In 
the  limiting  case  of  the  semi-infinite  needle  2a  -  0  the  motion  is  again 
of  the  diffusor  type.  A  comparison  with  the  analogous  result  presented 
in  [10]  reveals  that  the  velocities  of  plane  and  axisymmetric  diffusor  and 
jet  flows  with  radial  stream  lines  depend  on  the  polar  angle  9  in  exactly 
the  same  manner.  While  the  complete  Navier-Stokes  equations  yield  a  so¬ 
lution  which  represents  a  strictly  radial  diffusor  or  jet  flow  in  the 
plane  case  (see  the  Jeffery  -  Hamel  solution  in  [8  and  10]),  no  such 
integral  exists  in  the  axisymmetric  case. 

For  exact  solutions  of  the  Navier-Stokes  equations  of  the  form 
ijr  =  rF(9)  see  the  papers  of  H.  B.  Squire  [12  and  13]  . 
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TABLE  1:  Eigenvalues  \  =  ^  of  Ax i symmetric  Slow  Motions  Past  Conical  Surfaces  (12°  <,  la  <  160°) 
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